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Abstract 

Starting from known solutions of the functional Yang-Baxter equations, we exhibit 
Miura type of transformations leading to various known integrable quad equations. We 
then construct, from the same list of Yang-Baxter maps, a series of non- autonomous 
solvable recurrences of order two. 

1 Introduction 



A number of discrete integrable systems have been described in the literature in the last 
decades. They are mostly discrete versions of known integrable differential equations and 
partial differential equations, retaining fundamental features of integrability, as for example 
the existence of a sufficient number of invariant tori and hierarchies of associated integrable 
equations. Important classification results [1] have been obtained, both for the ordinary dif- 
ference equations case (in particular integrable maps [2] and discrete Painleve equations [3]) 
and for the partial difference equations (integrable lattice equations and Yang-Baxter maps 
among others [H EJ E]). 

We show that there exist bridges between different systems. We examine two such sit- 
uations. Starting from Yang-Baxter maps, given as lattice equations on a quad graph with 
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variables on the bonds, we first construct (section ED, via Miura type of transformations, 
quad maps with variables on the vertices, establishing a correspondence between the classi- 
fications of jl] and [5]. We then construct, from the same list of Yang-Baxter maps a series 
of solvable recurrences (section HD . 



2 The setting 

The underlying geometry behind the models is a two dimensional square lattice with variables 
associated to the edges and/or to the vertices. Each elementary cell is labelled by a pair of 
integers (m, n), which is the pair of coordinates of the lower left corner of the cell. Whenever 
possible we will use the standard notation "bar" and "tilde" to indicate shifts in the n 
and m direction respectively i.e. w = w(m,n), w = w(m,n + 1), w = w(m + l,n) and 
w= w(m + 1, n + 1). The elementary cell is depicted in Figure 1. 



x 



U 

Figure 1: Original cell 

The variables (x,y,u,v) are assigned to the bonds as shown in Figure 1. The models 
are defined by giving a 2 — > 2 map (x,y) — > (u,v), verifying the functional Yang-Baxter 
equations (aka Yang-Baxter maps). 

We start from the models described in [TJ E] and use the same labelling, namely the "F" 
and "H" families, which we reproduce here: 



u = ay P, v = (3x P, 



u = lp, v = ^-P, 
a 



u = -r, v 
a 



x 
J 



u = y P, v = x P, 



u = y + P, v = x + P, 



P 
P 
P 



(1 - P)x + - a + (a - l)y 
0(1 — a)x + (a — (3)yx + a(/3 — l)y' 
ax — f3y + (3 — a 



x-y 
ax — (3y 



x-y 

„ 8 — a 
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x-y' 

a — 8 
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(Fn) 
(Fin) 
(Fjv) 
(Fv) 
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= x Q, 
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u = 


y-P, 


V 


= X + P, 


P 



(l-(3)xy + ((3- a )y + f3(a-l) 
(1 — a)xy + (a — (3)x + a(f3 — 1) ' 
a + (/3 — a)y — (3xy 



f3 + (a — (3)x — axy' 
ax + (3y 



x + y 
axy + 1 
ftxy + 1 
a — [3 
x + y 



(Hi) 

(HA) 
(H^) 
(Hy) 



We will discard the cases HmA and Hy which are equivalent respectively to Fm and Fy, 
through the change u — >• —u and v —> —v. 

In section [3] we establish a correspondence between the original models and known inte- 
grable "quad equation", that is to say partial difference equations, which happen to fall in 
a classification given by the same authors [5]. 

In section |4] we construct from the same family of models a number of solvable ordinary 
difference equations. 



3 A correspondence via Miura transforms 

We will bring in vertex variables (w, w, w, w) (see Figure 2), subject to some local condition, 
via a Miura type of transformation of the original variables (x,y,u,v). 
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x " 
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Figure 2: Decorated cell 



3.1 Fi 

The first case we examine is Fi\ 

u = a y P, v = f3 x P, 



P 



(1 -f})x + P -a+ (a - l)y 
(3(1 — a)x + (a — f3)xy + a((3 — l)y 



(1) 
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Define vertex variables (w, w, w, w) via the Miura transformation: 

x = pww, u = pww, v = qww, y = qww (2) 

It is straightforward to eliminate the x, y, v, u variables using Eq. ([1]) and obtain an equation 
for w alone. Setting a = p 2 and (3 = q 2 , we obtain exactly equation A2 of [5]: 

(q 2 — p 2 )(wwww + 1) + p(l — q 2 )(ww + ww) — q{l — p 2 ){ww + ww) = (3) 

in canonical form. 

3.2 F n 

In the case of Fn 

u = -P, v = -P, P = (4) 

a p x — y 

we introduce the Miura transformation: 

x = pww, u = pww, v = qww, y = qww (5) 
with q 2 /p 2 = a/j3. We get 

~ p q 

q (ww + ww) — p (ww + ww) H = (6) 

q p 

which is equation H3 of [5]. 

3.3 Fni 

For system F m we start from 

V r> x ax - f3y 

u=-P, v = -P, P = (7) 

a p x — y 

and introduce the Miura transformation 

x = pww, u = pww, v = qww, y = qww (8) 

Substituting into the two equations we find the condition q 2 /p 2 = a/(3 whereupon w satisfies 
the equation 

q (ww + ww) — p (ww + ww) = (9) 

This is precisely the H3 equation of [S] for the special value 5 = 0. A more familiar inter- 
pretation of this equation is possible. Indeed solving for w we obtain 

qw-pw 

w = w — (10) 

pw — qw 

which is just the potential form of the discrete modified KdV equation. 
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3.4 Fiv 

The defining equation is 

a - B 

u = yP, v = xP, P = l (11) 

x-y 

Introduce the Miura transformations 

x = p + w + w, u = p + w + w, v = q + w + w, y = q + w + w (12) 

where B — a = 2(p — q). Eliminating the x,y,v,u variables using Eq. ( ITT]) we obtain an 
equation for w alone 

(w — w)(w — w) + (p — q)(w+ w + w + w) + p 2 — q 2 = (13) 

Equation (TIB"]) is precisely equation H2 of [5]. 

3.5 Fy 

The defining equation is 

a - B 

u = y + P, V = x + P, P= (14) 

x-y 

We introduce the w variable by 

x = w + w, u = w + w, v = w + w, y = w + w. (15) 

It is straightforward to eliminate the x, y, v, u variables using Eq. fll4p and obtain an equation 
for w alone. We find 

(w — w)(w — w) = a — 8 (16) 
which is the potential form of the discrete KdV equation that is to say equation HI of [5]. 

Next we turn to the study of the H family. 

3.6 #i 

The first case we examine is Hi 

u = yQ-i v = xQ g= (l-^ + (/3-a)y + /3(a-l) ( } 

The Miura transformation is 

x = pww, u = p{ww)~ 1 ) v = qww, y = q{ww)~ 1 (18) 

A direct substitution into Eq. ( 1TT1) leads to an equation which does not have the canonical 
A2 form. However introducing a = p 2 and 8 = q 2 we obtain 

(q 2 — p 2 ) (wwww+ 1) + p(l — q 2 ) (ww + ww) — q (1 — p 2 )(ww + ww) = (19) 

which is A2 in canonical form. 
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3.7 Hu 

The defining relation is 

n -i n n a + {fi - a)y - /3xy 

u = yQ , v = xQ, Q = — — — 20 

p + {a — p)x — axy 

The Miura transformation is 

x = ww, u = (ww) -1 , v = ww, y = (ww) -1 (21) 
and substituting into Eq. (12"01 we obtain 

(a — f3)(wwww + 1) + f3(ww + ww) — a(ww + ww) = (22) 

This equation could be interpreted as a special form of the A2 equation of [5]. Since 
Eq. ()3]) has the generic form of A2 it is convenient to start from this form and put q = 1 + ea, 
p = 1 + e/3. Taking the limit e — > we find precisely Eq. (122]) which establishes its relation 
to A2. 

3 . 8 -f^ni 

We start from 

V n n ax y + 1 fo>}\ 

w=7v v = xQ, Q = — (23) 

Q pxy + 1 

The Miura transformation is now 

x = pww, u = —p(ww)~ 1 , v = qww, y — — q (ww)~ (24) 

Eliminating x, y, v, u and introducing ap 2 = 1, (3q 2 = 1 we obtain the equation 

q (ww + ww) — p (ww + ww) = (25) 

which is equation H3 of [S] for the special value 5 = 0. 

We have thus established correspondences via appropriate Miura transformations, be- 
tween the two classifications of E] and [S] . 

Remark: All these constructions can be made non-autonomous [8l |9| [TO]. 



4 A second dynamical interpretation 

One may view the original lattice in a different way: the bond variables may be assigned to 
the vertices of a finer square lattice, as shown in Figure 3. 
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Figure 3: A change of point of view 



We get a new square lattice whose elementary square cell is, after rotation of — ? : 



I'- 



ll, 



Figure 4: New cell 



For all the models described in [7] the conditions relating u, v, x, y allow to define various 
rational 2x2 maps: they give u and v in terms of x and y, but also x and y in terms of 
u and v, as well as u and y in terms of v and x, or v and x in terms of u and ?/. This is 
the "quadrirationality" property. We obtain in such a way four 2x2 maps, defined for each 
elementary cell. These cells form a new two-dimensional lattice. Notice that there are two 
kinds of cells in this lattice. Half of them contained a vertex of the original lattice, and the 
other half did not. In the original lattice, we had relations only for the latter kind of cells. 
We will impose the existence of relations as above for all the cells, allowing non constant 
parameters. 

This will lead to a possibly non autonomous recurrence of order two. The variables, 
denoted by gned to the vertices of the new lattice. Since on each cell we have 

four (m, n) dependent 2x2 maps, we may, once given initial conditions £o,o and xo,i, calculate 
x over the whole lattice, provided some compatibility (zero monodromy around each vertex) 
condition is fulfilled. 
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Xm.n— 1 



■£777+1,77+1 



•^771,71—1 



771— 1,77— 1 



•^771—1, 



•^771— 1,71+1 



Figure 5: Local compatibility condition around vertex (m, n) 

We use the coordinates of the lower left corner to label the cells. Starting from x m>n and 
x mj n-i, we compute £ m ,n-i and x m+ i tn using the same F and H equations as in the previous 
section, with (m,n) dependent set of parameters r mj „. Going around the vertex (m,n) we 
calculate the surrounding values of x. This yields after four steps a value for x m ri _i which 
has to coincide with the initial condition. 

Remark: Doing so we reduce the dimensionality of the system, going from partial dif- 
ference equation to ordinary difference equation. The reason is that we impose relations 
between variables which were not related in the original model, actually exactly on half of 
the cells of the new lattice. 

In what follows, we shall show that it is indeed possible to satisfy the compatibility 
constraint. Moreover the resulting equations are explicitly solvable ones and given that the 
compatibility constraints are also explicitly solvable this allows us to give the solution for all 
x's in closed form. 



4.1 Fi 



In the case of Fj, the starting equations are 



u = ay P, v = (3 P, P 
which we rewrite as 

•^777+1,71+1 < -^777,7l2'777,n+l 



(1 - f$)x + P - a + (a - l)y 



/9(1 — a)x + (a — 0)xy + a(/3 — l)y 



771+1, 7 



/^m,n(l (^m,n)X m n -\- {ot m n ftm,n)X m ^ n X mtn -\-\ -\- ttm,n (/3t77,77 \^X m ^ n -\.\ 



@m,n(X (^777,77 )Xui, n ~\~ \PLm,n /3m,n)%m,nXm,n+l 0^m,n(^Pm,n ^-)X m ^ n 



+ 1 
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The compatibility conditions for the a and /3 are 

(26) 

and 

QWfl,n = firn,nfirn,n+lR, Pm+l,n+l = <^m.,n a; m,n+l-R (27) 

where 

^ _ («m,n ~ l)(Qm,"+l ~ 1) ~ (Pm,n ~ l)(^m,n+l ~ 1) 

The set of compatibility equations is itself overdetermined but it turns out that it is 
consistent and the solution is periodic in both m and n. We have 

«3m,3n = «3m+l,3n+l = «3m+2,3n+2 = OLq 

«3m,3n+l = «3m+l,3n+2 = «3m+2,3n = OL\ 
«3m,3n+2 = «3m+l,3n = «3m+2,3n+l = a 2 
03m,3n = 03m+2,3n+l = /^3m+l,3n+2 = A) 
03m,3n+l = P3m+2,3n+2 = P3m+l,3n = /?1 
03m,3n+2 = /33m+2,3n = /33m+l,3n+l = 02 

where a 2 , /3 2 are given in terms of a , cm, /3o, /?i using Eq. (T27j). The solution for x m>n is 
equally periodic, with period 3 in both m and n directions. In this case however all nine 
X3m+j,3m+k for j, k — 0,1,2 are generically different. 



4.2 F^ 



II 

We start from the equations 

V r> x td td ax - f3y + (3 - a 

u = -P, v = -P, P = (29) 

a p x — y 

which we rewrite, introducing 7 = a/ (3, as 



•£m+l,n+l •Em,n+1 



The compatibility condition is now 

^m+l,n+l^rn+l,n'~1m,n-\-l'yrn,n 7«2+l,n7«'i,«+l'")Vrt,n r )m-\-l,n-\-l r )m-\-l,n'~1m,n+l ~l - 7«}+l,^ ~^~ n fm,n+l 1 

which can be rewritten as 

(1 - 7 m+ i >n+ i)(l - 7 m>n ) = (1 - l/7 m+ i >n )(l - l/7 m>n+ i) (31) 
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which can be recognized as the linearisable equation of Hydon-Viallet (Eq (22) in [IT], 
originally given in terms of the variable z — 1 — 7)- The solution to Eq. (1311) can be obtained 
most easily if we introduce an auxiliary variable x an d express 7 as 

(1 Xm+l,ra+l)(l Xm,n) f r i r )\ 

7m,n — JZ \7i \ \ 6Z, > 
U — Xm+l,n){± — Xm,n+l) 

with x satisfying the equation 

Xm-\-l,n+lXm,n Xm - 1 .nXrn.n — 1 

(33) 

i.e. x is t ne product of two free functions Xm,n = f( m )g( n )- Using this solution we can 
obtain the general solution for 

_ af(m)g(n) + bf(m) + cgjn) + a + 1 
1 - f{m)g{n) 

where a, b, c are constants related by a(a + 1) = be. 



4.3 F m 

In the case of Fm, where 

u = -P, v = -P, P = (35) 

a p x — y 

we introduce the parameter 7 = a//3, and rewrite the equations as 

(•£771,71 ^m,n+l/7m,n 1 _ f 1m,n%m,n %m,n+l \ (na\ 

•£771,71 •£m,n+l / \ •^m,n ■^m,n+l / 

The compatibility condition is exactly the same as before, namely 7 satisfies equation (|3ip . 
In this case, it is more convenient to express 7 in terms of an auxiliary variable u as 

^771+1,71+1^771,71 /orr\ 

7m,n = (37) 

^771+1,71^771,71+1 

and take u satisfying the equation 

w m+l,n+l + w m,n = w m+l,n + w m,n+l 

(38) 

i.e. u is the sum of two free functions u m ^ n = f(m) + g(n) which leads to the solution for 7 

(/(m + l)+g(n+l))(/(m)+g(n)) 
7m ' n (f( m +l)+g(n))(f(m) + g(n+l)) { ] 

The general solution for x can now be given 

_ (f(m) + b)(g(n)-b) 
m ' n ~ f{m)+g{n) m 

where a, b are constants. 
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4.4 Fiv 



The defining equations 



u = y P. v = xP, P = l- - — - (41) 

x-y 



are rewritten as: 



/ -i lfm,n \ / 7m,n \ ( AO\ 

•^m+1,71+1 %m,n+l I I > -£771+1,71 %m,n I J- I l^^J 

\ •Em,n -£y7i,n+l / \ -£771,71 -£m,n+l J 

where 7 = a — /3. 

The compatibility condition in this case turns out to be 

r )m+l,n+l r Ym,n 7771+1,7i7t71,71+1 

(43) 

the solution of which is 7 mjn = <j)(m)ip(n), where and ip are free functions of their argument. 
At this point is is convenient to introduce two functions f(m) and g(n) such that f(m + 
1) — f(m) = <f)(m) and g(n + 1) — g(n) = ip(n). The general solution here can be written as 

x m ,n = (/(m) + a)(g(n) + b) (44) 

where a, b are constants. 

4.5 F v 

The defining equations are 

u = y + P, V = x + P, P = (45) 

x-y 

with the same definition of 7 as before, and we rewrite them as 

, Im.n 7m,n (A(\\ 

•^m+1,71+1 -£771,71+1 I ) 2-771+1,71 %m,n ~T V^'-'J 

•^m,n • L "m,n-\-l Jj m,n J -"tn,n+l 

The compatibility condition in this case turns out to be the same as in the previous case, 
namely Eq. (143]) . Using the same functions / and g as above, it is straightforward to verify 
that 

x m ,n = af(m) + ^ + b (47) 
a 

is the general solution of system (1461) . where a, b are constants. 
Next we turn to the study of the H family. 
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4.6 Hi 



The first case to examine is Hi 

u = yQ^ v = xQ Q= (l-^ + (/3-a)y + /3(a-l) ( } 

which we rewrite, after a translation a — 7 + 1, /3 = 5 + 1, as 

Om,n-Em,n-Em,n+l ~T ^7m,n " m,n )-^m,n+l '~1m,n\'Jm,n < J- J 

_ ^m, n^m,n^m, n+ 1 "I - (7m,n ^m,ri) -^771,71+1 Tm,n('^m,ri "I - 1) fcrA 

%rn+l,n = %m,n 7 7 7 7 7 \ l^^J 

7111,11^111,^^1,11+1 ^7m,n "m,,nJ3'm,n "m,nl7 m , n ~r J- J 

In this case the two parameters 7 and 5 cannot be combined into a single parameter. The 
compatibility condition turns out to be 

7m+l,n $m,m ^m+l,n 7m,n (^1) 

which means that 7 and 5 can be expressed in terms of two free functions /, g of n alone: 
l2m,n = S 2m +i, n = f(n), 72m+i,n = hm,n = gin). The solution of Eqs. (j49|), (EDJ) is also 
periodic in the m-direction: x m+ 2,n = Xm,n- Thus the equation becomes a one- dimensional 
one for the variable x at given parity of m. Once we have obtained the solution for one parity 
the solution at the other parity is obtained through Eq. (149]) or (150]) . Taking the index m to 
be even we obtain the equation for x 

(g(n)g(n - 1) - f(n)f(n - l))x n _ix n x n+1 + 
+ g(n)(f(n - 1) - g(n - 1)) x n x n+1 + f(n - l)(f(n) - flf(n)) x„x„-i + 
+ /(n - l)(<?(n - 1) + l)(/(n) - g(n)) x n+1 + g(n)(f(n) + l)(/(n - 1) - <?(n - 1)) + 
+ - l)(/(n) + /(n - 1) + 1) - f(n)f(n - l)(g{n) +g(n-l) + l) = (52) 

where we have dropped the even index m. (In fact the equation for odd-m is the same up 
to the permutation of / and g). Equation (l52~l) is linearisable. A study of the growth of the 
degrees of the iterates on an initial condition, using algebraic entropy techniques, leads to 
the following sequence of degrees: 0,1,1,2,2,3,3,4,4,5, .... This suggest that the mapping 
defined by (152]) is of Gambier type [12] and this turns out to be indeed the case. We find 
the invariant 

K _ , , )n f( n )( x n + l)(zn+i - 1) + g{n)(x n - l)(z n +i + 1) -2f(n)g(n) , , 
1 } (f(n)-g(n))(x n -l)(x n+1 -l) 1 1 

which allows to reduce Eq. (15"2]) to a homographic mapping. 



4.7 Hu 

The defining equations 



u = yQ _1 , » = »0, (54) 

p + {a — p)x — axy 
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can be cast into the form 

1 \ fm,n l)-£m,n '~Ym,n-Em,n-Em,n-\-l fz.z.\ 

%m+l,n+l = x m,n+l T"7i \ V / 

7m,n i V 7tt7,77,)'£777,77+1 •^m,n-^m,n+l 

_ r Ym,n (1 Om,n.)'£m,n+l 2'777,77'2'777,77+l 

•^rn+1,77 •Em,n~ j j 7 77 v."-*"/ 

where 7 = a//3. It turns out that the compatibility condition for 7 is precisely Eq. fl3Tj) i.e. 
the Hydon-Viallet equation. We give the solution for 7 in the same form as in the case of 
Fji, namely 

(1- f(n + l)g(m + !))(! -f(n)g(m)) 
^ (l-/( n +l) s (m))(l-/( n ) ff ( m +l)) 1 1 

We can now give the solution for x which turns out to possess an even-odd freedom 

x = (1 + a n f(n))(g(m) + a n+1 ) 
(1 + a n+1 f(n))(g(m) + a n ) 

where a n is a period-two constant, i.e. a n+ 2 = a n . 



4.8 



The set of equations H 



111 



V n n ax y + 1 ftia\ 
u = — , v = xQ, Q = — (59) 

Q (3xy + 1 



yield 



•£777+1,77+1 *£t77,77+1 j 7? %m+l,n %m,n~^ j 7" 

The latter can be written in more compact form if we introduce a = 7 + 5 and /3 = 7 — 5. 
We obtain 

^777, 77^777+1, 77+1 ^771, 71 + 1^777+1, 71 (^l) 

(7777+1,77 _ 7777,rt)^m + l,77+l + (7777 + 1,77+1 — 7777,77+l)^777 + l,77 = (62) 

which can be readily integrated to 

CW = f(n){h{n) + (-l) n k(m) + g{m)) (63) 

f3 m ,n = f{n)(h(n) + (-l) n k{m) - g{m)) (64) 
Next we remark that from Eq. (|60|) we have 

^777, 77^ 777, 77+1 = ^777+1,77 ^777+1,77 + 1 = 4 > ( n ) (65) 

where <fi is, obviously, independent of m. Given the expressions of a and f3 we can obtain the 
value of 1/0 = f(n)(—h(n) + ( — l)"c), where c is an integration constant, the value of which 
is fixed by the initial conditions. The values of x can now be obtained from the recursion 

_ 0(n + l) 

•^777,77+2 1/ \ -^777,77 (00 1 

<j>{n) 
13 



Using the expression for we can, starting from some initial condition x ,0) ^o,i obtain all 
£ ,n for all n. Once the x 0:n are obtained we can compute x m>n by 

_ c-k(m) -g(m) 
c — k{m) + gym) 

c - kjjn) + gjjn) 

x m+l,2n+l — T7 \ ^ x m,2n+l {vo ) 

c — kym) — gym) 

Thus the solution to Eq. (160]) can be expressed as a double semi-infinite product. 



5 Conclusion 

Integrable discrete equations appear in an innumerable number of avatars, and we have 
exhibited some relations between such forms. 

The relations we have described enter the general problem of classifying integrable dis- 
crete equations. One should nevertheless emphasize that the correspondences we have es- 
tablished are not invertible. They consequently do not define equivalences between models. 

This is reminiscent of what happens when one goes from a discrete equation to its "po- 
tential form", or when one performs reductions by symmetries. These are other types of 
transformations establishing relations between different integrable systems, and the links 
between these various correspondences and the one we presented here still need to be ex- 
plored. 

Acknowledgment: We would like to thank Jarmo Hietarinta for stimulating discussions 
in the early stage of this work. 
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